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Abstract. The electronic structure (electron charges and spins) of a per- 
fect crystal under external magnetic field is analyzed using the spin-polarized 
Thomas-Fermi-Dirac-von Weizsacker model. An extension of the classical 
Cauchy-Born rule for crystal lattices is established for the electronic struc- 
ture under sharp stability conditions on charge density wave and spin density 
wave. A Landau-Lifshitz type micromagnetic energy functional is derived. 



1. Introduction 

The present paper is the fourth of a series of papers that are devoted to the 
study of the electronic structure of smoothly deformed crystals or crystals in an 
external field, by analyzing various quantum mechanics models at different levels 
of complexity, including the Kohn-Sham density functional theory, Thomas-Fermi 
type of models and tight-binding models. Our overall objective is to establish the 
microscopic foundation of the continuum theories of solids, such as the nonlinear 
elasticity theory and the theory of magnetic materials, in terms of quantum me- 
chanics and to examine the boundary where the continuum theories break down. 

In previous work [5-7], we have studied the nonlinear tight-binding model and 
the Kohn-Sham density functional theory, and established their continuum limits 
for smoothly deformed crystals. As a byproduct, we also derived macroscopic mod- 
els for the piezoelectric effect of a material from many-body quantum theory. In the 
present work, we will focus on the magnetic properties of a material and the associ- 
ated spin waves. We choose the spin-polarized Thomas-Fermi-Dirac-von Weizsacker 
(TFDW) model as our starting point. This is a representative of a simplest version 
of the density functional theory in which the electronic structure is represented 
solely by the electron density instead of the wave functions as in the Kohn-Sham 
theory. As in the previous work, our reference point is the electronic structure of 
a perfect crystal, i.e. the equilibrium crystal lattice without external fields. The 
solution of the TFDW model we are interested in is a continuation of the solution 



Date: June 17, 2010. 

The work of W. E was supported in part by the NSF grant DMS-0708026, grant DMS-0914336, 
the ONR grant N00014-01- 1-0674 and the DOE grant DE-FG02-03ER25587. J. Lu would Uke to 
thank Professor Robert V. Kohn for helpful discussions. 

1 



2 



WEINAN E AND JIANFENG LU 



for the unperturbed system. This works as long as certain stability conditions are 
satisfied. One main objective is to identify these stabiUty conditions, with emphasis 
on the stability of spin density waves (magnons). The overall strategy is similar to 
the one in the previous paper [6] and [5]. 

Another related work that we should mention is [1] where the Thomas-Fermi- 
von Weizsacker model (without spin-polarization) was studied for smoothly de- 
formed crystals. It was shown that in the continuum limit, the total energy of the 
system converges to a limiting value given by the extended Gauchy-Born construc- 
tion. There, the stability condition is automatically satisfied since the model is 
convex. [4] extended this kind of results to the tight-binding models. The strategy 
in the current series of papers is quite different from those in [1] or [4]. Here our 
emphasis is on the stability conditions. 

Throughout this paper, wc use the notation < for inequalities up to an absolute 
constant: / < g if / < Cg where C is an absolute constant. Sometimes, it is more 
convenient to explicitly use C to denote the constant, which might change from 
line to line. When it is necessary to specify the dependence of the constant on 
parameters, we will use the notation C(a, b) to indicate that the constant depends 
on parameters a and h. 

Standard notations are used for function spaces like L'p , and W^''^ so on. We 
also need function spaces for periodic functions. Let L C be a lattice with unit 
cell r. Denote the reciprocal lattice as L* and its unit cell (the first Brillouin zone) 
as r*. For a given n e N, we define 

(1.1) LP = {/ e ^'(R3) I Tfl/ = /, Vi? G nL; / da; < oo}, 

JnT 

with norm given by 

(1-2) \\f\\i^.^={n-^ jjn^dx 

and similarly for L^. Here tr is the translational operator with translation vector 
R-i ijRf){x) = f{x — R). The periodic Sobolev space is defined similarly 

(1.3) H^ = {f& I TRf = f,\/R€ nL; / e H>'{nT)}, k€Z+ 

with norm 

(1-4) ii/iks = E w^yhi, 

\a\<k 

where a denotes a multi-index and the corresponding partial derivative. More- 
over, we also define the periodic homogeneous Sobolev space with index —1 as 

(1.5) H-\nT) = {fe^'{R')\TRf = fyRenh, ^ J_|/(fc)|2 < oo}. 
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Here, {/(fc)} denotes the Fourier coefBcients of the nF-periodic function / 
(1.6) 7(/e) = (27r)-3/2 / f{x)e-"'--dx. 

JnT 

The space H~^{nT) is a Hilbert space with inner product 
(1-7) (/.5)H-i(„r)=47r jLWm- 

For Banach spaces X and F, ^{X,Y) denotes the class of bounded hnear op- 
erators from X toY and |H|^(x,y) denotes the operator norm. When X = F, we 

Vise^{X)=^{X,X). 

2. Spin-polarized Thomas-Fermi-Dirac-von Weizsacker model 

We consider the spin-polarized Thomas-Fermi-Dirac-von Weizsacker (TFDW) 
model. We will restrict ourselves to the collinear case, that is, we assume that the 
applied magnetic field is parallel to a fixed axis (with possibly varying amplitude). 
In the spin-polarized TFDW model, the electronic structure is characterized by 
the spin-up and spin-down densities, corresponding to the density of spin-up and 
spin-down electrons. Denote by p+ and p_ the spin-up and spin-down densities 
respectively, the energy of the system is given by 

(2.1) i\p^,p-) = j pT + p'-1' + \^Vp^? + \^^?- j pT + p1' 

+ \^^P - Pb^P- Pb) - j hm. 

Here p = p+ -I- p_ is the total charge density, m = — p_ is the spin density 
or magnetization density, and pf, is charge contribution from the nuclei and core 
electrons. The shorthand notation D{-, ■) is defined as 



\x - y\ 

Finally, h is the external magnetic field, which is a scalar, since we have assumed 
the coUinearity. 

The density satisfies the normalization constraint 



(2.2) J p = J p,, 

and the positivity constraint p+, p- > 0. To deal with the positivity constraint, it 
is often convenient to introduce the variables 

^+ = Vp+, ^- = sfP^- 
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In terms of and v- , we have 



(2.3) u.) = I ,.1'/' + v'y^ + |Vi.+ |^ + |Vi.-|^ - j v'f + 



+ ^D{p- pb,p- pb)- I hm, 



where 



2 I 2 

p = u_^_ + u_, 



The energy is given by the variational problem 

(2.4) E'' = inf I^(v+,v_) 

with the normalization constraint (2.2). The Euler-Lagrange equation associated 
with the energy functional reads, 

(2.5) - Ai/+ + - + Vv+ - hv+ = Az^+; 

(2.6) - Az/_ + §1/!./^ - |zy!./^ + Vu_ + hv_ = 

(2.7) -Ay = 47r(p-p6), 

where A is a Lagrange multiplier corresponds to the normalization constraint. 

In this work, wc will consider the TFDW model for crystals with and without 
the external magnetic field. Wc denote L the underlying Bravais lattice of the 
crystal and F the unit cell of the crystal. The charge background pb is F-periodic 
and assumed to be smooth (in other words, we are taking a pseudo-potential ap- 
proximation). When there is no external magnetic field, the electronic structure is 
characterized by the periodic-TFDW model, with energy functional 

(2.8) /°(.+ , V.) = vf' + + |V.+ |^ + IV.- 1^ - .'f + 



+ \dt{p- Pb,p- Pb), 



where the periodic Coulomb interaction is given by 

Drif,g) = (.f,5>H-i(r)- 
The electronic structure is given by the variational problem 
(2.9) M 

with the normalization constraint that 



T 

Here Z is some fixed integer. 
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The Euler-Lagrange equations associated with (2.9) arc given by 

(2.10) - Au+ + f i.^/' - 14/' + Vu+ = 0; 

(2.11) - Ai/_ + fz/I/-' - |z/!/'* + Vu_ = 0; 

(2.12) - Ay = 47r(p-pb). 

The equations are defined on F with periodic boundary conditions. Note that 
we have absorbed the Lagrange multipher into the potential V, as (2.12) only 
determines F up to a constant. The normalization constraint is also implicitly 
imposed as the solvability condition of (2.12). 

It is easy to see that the minimum of the variational problem (2.9) is achieved 
and the minimizers satisfy the Euler-Lagrange equations. The minimizers might 
not be unique since the Dirac term (— / v^^^ + 1/!^) is concave. Let us take one of 
the minimizers, denoted as z^+,per and t'_,per) and let us denote the corresponding 
potential as V^er (with the Lagrange multiplier included). By standard elliptic 
regularity theory, it is easy to see that J^±,por, ^pcr <= C°°(r). We extend them 
to the whole periodically. It is also straightforward to see that v±^-pcr are non- 
negative. It is possible that z^_^per = is a minimizer (while the corresponding 
i/+,per > 0). To exclude such cases, we will assume that there exists a positive 
constant such that 

(2.13) i^±,per(a;) > C,, x& T. 

We may also consider the energy functional defined on the supercell nV for n e N: 

(2.14) IIAV+,U_)=[ ^|°/3 + ^iO/3 + |V^+|2 + |V,._|2- / i.«/3 + ^8/3 

JnV JnT 

+ ^Dnr{p- Pb,P- Pb), 

where 

D„r{f,9) = {f,9)H-^nry 
The corresponding Euler-Lagrange equations have the same form as (2.10)-(2.12) 

defined now on nT with periodic boundary conditions. The functions i^i^per, Vper 
(recall that we have extended them to periodically) still satisfy the Euler- 
Lagrange equation and hence are stationary points of the functional (2.14). 

3. Stability analysis 

3.1. Stability of Electronic Structure. For any n € N, let us define the linear 
operator C, which is the linearization of the Euler-Lagrange equations (2.10)-(2.12) 
on the domain nT, given by 



(3.1) 



UJ+\ 













H 





£_ 




wj 




w 
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where the operators £± are given by 

(3.2) L±bJ = —Aw + ^i'±^u) — -^I'j^^co + Vu. 



Proposition 3.1. Assume y±, V G W^'°°, for some k gN. Then C is a bounded 
operator from (iJ^+^)^ to {H^)^, 

Moreover, if v±, V G , then £ is a self-adjoint operator on the domain 'D{C) = 
Proof. Denote 



then we have 



. A , 4/3 20 2/3 , , „r 

/_ = - Ac.. + ^v^J''w. - ^uTuJ- + Vuj. + Wu_; 

It follows that 

< C(||i/+||^;^fe,oo, ||y||^;^fe,oo)||a;+||^fc+2 + ||!/+||^fc,oo||T4^||jjfe+2; 

WqWh^ < ^W\\h'^+'2 + ||z/+||,^fc,oo||w+||^fc+2 + ||z/_||^fc,oo||a;_||^fc+2. 

Hence, £ is bounded with the desired estimate on the operator norm. 

The self-adjointness of C is an easy consequence of the Kato-Rellich theorem [9] , 
since v± and V viewed as multiplicative operators on H'^ are infinitesimally small 
with respect to the Laplacian operator. □ 

Consider specifically the case when u± and V given by the unperturbed system: 







\ '^+,per 





J^+,per\ 


(3.3) 


c -1 

■*-'per — 1 





r 

— ,per 


^— ,per 1 






\^'+,per 


^— ,per 





with £±,per defined by 

(o A\ r A , 4/3 20 2/3 , X, 

(3.4) £±,perW = -Aw + Yi^±',per'^ " y'^i.per'^ + ^perW- 
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Assumption A (Stability of ttie electronic structure). There exists a constant M 
independent of n, such that for any n, 

ll'^pcrl|j5f((L2)3) < M, 

or equivalently, since £per is self-adjoint 

dist(0,spec(£per)) > 1/M. 

Under the stability assumption, we can actually obtain estimates of C~l^ acting 
on higher order Sobolev spaces. The following result is standard, we include the 
proof here for completeness. 

Proposition 3.2. Under Assumption A and assume u±, V G W^'°°, for some 
k gN. Then we have 

ll'^pirll^((H^)3,(H^+2)3) < C{k)M. 

Proof. Let us consider A; = first. It suffices to prove the estimate 

(3.5) P£->||(ff2)3 <M||w||(i.)3, 

for any u = {uj+,uj-,W) € (-L^)^, where A is the operator 

-A 

A=\ -A 



The left hand side of (3.5) equals to 
/-A 

-A 



where 










£pg> = u- 









\^-|-,per 


^— ,per 



-,per I *-per'- 



P^_ 35 4/3 _ 20 2/3 y 
^± — -9"^±,per -9"^±,per + ''per- 

Therefore, 

\\A^V^M\{LIY ^ +max(||F±||Lc»,j|z/±,poi.||L-)||/:p'erl|.5f((L=)3)l|w||(L=)3 

<M||m||(l2)3. 

Suppose the statement of the Proposition is proved for k < ko, let us consider 
k = ko + 1. Since 

^ -^per ■'^per * ^ L ^ ' -^perj 

— r~^T7 — r~^\T7 c 

'-'per " '-perl * ' '-'perj-^per) 

it suffices to control 



ll'^peit^' ■^Per]'^perll^((JJ*)3_(H^+i)3), 



8 



WEINAN E AND JIANFENG LU 



where V = /3 V with I3 the 3x3 identity matrix. Note that 

ll'^per[^5-^per]'^pe^rll.Sf((Hfe)3^(H^+i)3) < \\^per\\^{{H^-^)3 ^(H^+^)3) 

^ II ■^per] ll^((if^+i)3_(Hfe-i)3) ll'^perll.5f((ff^~^)3,(JJ^+i)3)- 

By assumption, the Proposition holds for A: — 1 = fco, hence it suffices to control 
the commutator [V,£per]- An explicit calculation yields 

/ VF+ V//+,pe/ 



[V,>Cper] 



VF_ Vj^-,per I , 



\Vz^+,per Vl^_,per 

and hence the bounds follow from the regularity assumptions on t'i^per and Vper- 
The proposition is proved. □ 

The Assumption A is stated in terms of the operator £per acting on a series of 
spaces (i^)^. Using the Bloch-Floquet decomposition (see [10] or [7] for an intro- 
duction), we may obtain an equivalent characterization of the stability assumption. 
Note that £pcr commutes with the translational operator with respect to the lattice 
L since i'±,pei and Vper are F-periodic. Denote T* as the unit coll of the reciprocal 
lattice of L, the Bloch-Floquet decomposition of £per is given by 

(3.6) £per = f £j,per d^. 

Here for any ^ € F*, -Cj.per is the operator 

(3-7) ■^^.per = £__{_per 

\ ^+,per ^— ,per 

with £±,£,per given by 

(3.8) £±,perW = -AgW + ^U^l^^U - ^u'^^l^^U + FperW- 

We also have for any ^, the operator £{,per defined on the space (i|)^, where 
= {/ e LL(M3) I e'i-f{x) F-periodic}, 

is self-adjoint [10]. 

Using Bloch-Floquet decomposition of £per, the stability assumption can also be 
formulated as 

Assumption A'. There exists a constant M, such that for each ^ e F*, 

WClU^r^D^ < M. 

The proof of the equivalence of Assumption A and Assumption A' is parallel 
to the corresponding results in [5,6] and is standard from Bloch-Floquet theory; 
hence, we omit it here. 
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3.2. Example of stability and instability. Let us consider the jellium model as 
an example to understand better the stability assumptions. For the jellium model, 
the charge background is a constant function ph{x) = po- Define uq = \py^ , a 
solution to the Euler-Lagrange equations are given by 

v+{x) = vo, v-{x) = VQ, V{x)=Q, X = ^Vq'^ - ^vl'^. 

Therefore, 

/-Q n^ r r \ ^ '^^ 4/3 8 2/3 

(3.9) L+w = L^ijJ = —Aw + ~ g^^o 

We use Fourier transform to analyze the operator C 

(e + f4^' - 

(3-10) C,= \ ^2 + 1,^4/3 _ 8^2/3 

One eigenvalue of the matrix is 

(o ^^\ \ (^2 , 20 4/3 8 2/3 

(3-11) Aj,i = C + —1^0 - gZ^o > 

corresponds to the eigenvector (1,-1, 0)""^. The other two eigenvalues are given by 

r-1 2 20 4/3 8 2/ 



^'^19^ X 1 ^ 87r-1 ^2 , 20 4/3 8 2/3 

(3.12) A4,± = - — — ^ ^—v^ --V,' 



± 



\7^TV~~^2^^ 8 2/3\2 . 2 



Observe that these two eigenvalues correspond to action of on the subspace 
orthogonal to the vector (1,-1,0)^, and hence, the first two components are the 
same. 

It is easy to see that A^^i becomes positive if ^ is sufHciently large. To prevent 
A^,i from changing sign when ^ is small, we need 

fi , 20 4/3 8 2/3 ^ ri 

or equivalently 

(3.13) -0 > (|)'/'. 

For the other two eigenvalues, we have 

\ A 1 c2/t2 , 20 4/3 8 2/3\ 2 

The product becomes negative when ^ is sufficiently large, and hence we have 
Aj,_ < and Aj,+ > for ^ sufficiently large. To make sure they are nonzero for 
every ^, we need 

1 ^2/^2 , 20 4/3 8 2/3\ 2 w,- 
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which is equivalent to the condition that 

(3.14) - lu^/' > -8v^.o. 

This condition is weaker than (3.13). 

Hence, (3.13) guarantees that the three eigenvalues do not change sign for all ^, 
and hence the matrix is non-singular. 

Let us remark that physically, the condition (3.13) corresponds to the stability 
of spin-density-wavc, since Afe.i corresponds to the eigenvector (1,-1,0)"^, which 
increases (or decreases) the spin-up component, while decreases (or increases) the 
spin-down component, hence creates a spin-density- wave. On the other hand, the 
condition (3.14) corresponds to the stability of charge-density- wave, since the spin- 
up and spin-down components change together with the same amplitude. For the 
jellium case, since the condition (3.14) is implied by the condition (3.13), we observe 
that the spin-density- wave loses stability earlier than the charge- density- wave when 
the uniform background charge density is decreased. 

4. Cauchy-Born rule 

Let us first consider the case when the applied magnetic field is constant h{x) = 
h, and consider the cell problem 

(4.1) 7p\^+, u_) = J^ 4'/' + u'y' + + |2 - 1^ u'f + u'J' 

+ ^^r(/9 - Pb,p- Pb) - hmtot, 
where the periodic Coulomb interaction Dr is given by 

Dr{f,g) = (/,5)ij-i(r)> 
and the total magnetization mtot is given by 



mtot = 




Consider 

(4.2) E{h)= inf J^, 

with the normalization constraint j^p = Z. The Euler-Lagrange equations are 
given by 

(4.3) - A!/+ + Iv'f - + {V- h)v+ = 0; 

(4.4) - Ai/_ + IvJ^ - + {V + h)v_ = 0; 

(4.5) - Ay = 47r(p-p6); 

in the unit cell F, where his & constant. Let us denote the equations as 

T{u, h) = 0, 
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where u stands for the triple u = {v+, v-,V). 

The following theorem shows that the solution u exists provided the stability 
condition is satisfied and the constant applied magnetic field h is not too large. 

Theorem 1. Under Assumption A, there exist positive constants ho and S, and 
a unique C°° map from [—ho, ho] (Hf)^: h u{-;h), such that \\u{-;h) — 
Uper 11(^3)3 < S and 

J^{u{-,h),h) =0. 

Proof. We use the implicit function theorem. Let Uper be the triple (!^±,per, V^er), 
we have 

J-(Wper,0) =0. 

Let 6i be a positive constant to be fixed, consider the neighborhood around ttper^ 

T>= {u\\\u- Uper||(//3)3 < Si}. 

Denote u = {v±, V); by Sobolev inequality, we have 

Take sufficiently small such that \\v± — i^±,per||Lj° < Hence, for u G V, 

we have i'± > C„/2 > 0. It is then easy to see that viewed as an operator from 
2? X M3 {Hl f, J" is C°°. Notice that 

dJ-{u, h) 
6u 

has a bounded inverse from {Hl)^ to (-fff)^ due to the stability assumption and 
Proposition 3.2. Applying the implicit function theorem on J^, we arrive at the 
desired result. 

□ 

Remark. It is clear from the proof that we can consider solutions in spaces with 
higher regularity. The space (-ffi)^ is chosen here for proving the main result in 
the next section. 

We denote the solutions given by Theorem 1 as v^^cBi'', h) and i^-,cb('; h) respec- 
tively for the spin-up and spin-down components, and VcB{-]h) for the potential. 
Here /i is a parameter and z^±,cb and Vcb are F-periodic. 

For h & [—ho, ho], let Ch be the linearized operator around ucb(-, h): 

6T{u, h) 



= r 

U=UpeT,h=0 



(4.6) U . 

ou 



u=ucB{-',h),h=h 



given by 











(4.7) 

















Ch,+ i/+,cb(-;^) 

Ch,- i/_,cb(-;^) 

i'+,CB{-,h) u-,cB{-,h) 
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where the operators Ch,+ and Ch.- arc defined as 

35 20 

(4.8) Ch,+^ = -Aw + -^v+^^ix; h)ui - y'^+%b(.^' + (^CBix; h) - h)ui; 

(4.9) Ch-co = - Aw + yz^1(cb(^; f^)^ ~ y^-(cb(^; + {Vcb{x; h) + h)Lj. 

By Proposition 3.1 and using similar arguments as the proof of Proposition 3.2, we 
have 

Proposition 4.1. Under the same assumptions of Theorem 1, Ch as an operator 
from {Hi)^ to {Hl)^ is invertible, and the norm of the inverse operator is uniformly 
bounded for h G [—ho, ho] ■ 

Let us remark that Theorem 1 gives a map from h to the electronic structure 

(i^_i-,i^_) for the case when the external magnetic field is homogeneous. This is 
slightly difl:erent from the usual Cauehy-Born rule for crystals under deformation, 
where the strain is fixed - here h, the analog of stress, is fixed. One may consider 
the dual problem given by 

Ecnim) = inf ( / + i.'^' + |Vi.+ p + - / ul^' + 

v±:mtot=rn Jr 

+ ^^r(p-pfc,p-pb)|, 

where mtot = Jr'^+ ~ ^- constrained to be equal to m. This can be viewed as a 
Legendre transform of 

EcBih) = inf -f uf' + u'^' + |Vi.+ p + - j- + u'J'' 

\ , ,1 

+ -j^^^TKP - Pb, P - Pb) - i^/iTOtot. 

The formulation in terms of the magnetization m may bear more similarity with 
the conventional Cauchy-Born rule for lattices. 

5. Main results 

We turn to the situation when the system is under (a macroscopically heteroge- 
neous) external applied magnetic field. We will study the case when the applied 
potential is macroscopically smooth. The ratio of the lattice constant and the char- 
acteristic length of h will serve as a small parameter e. Given a fixed F-periodic 
function /i(-), two equivalent choices of scalings are possible: For any e = 1/n a re- 
ciprocal of positive integer, we may study a perfect crystal with applied field h{ex) 
in nV with periodic boundary condition; equivalently, we may rescale the system, 
so that the lattice constant becomes e and study the rescaled system with applied 
field h{x) in F with periodic boundary condition. We call the former choice the 
atomic unit scaling, and the latter choice the £-scaling. We will use atomic unit 
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scaling lor most part ol the paper, however, e-scaling is more convenient and is 
used for the two-scale analysis in Section 6.1. 

Under the influence of the external fleld, the electronic structure of the system 
is determined by minimizing the energy functional. 

(5.1) />+,._)=/ .W3 + ,iO/3 + |V.+ |2 + |V._|2_ / ,8/3^^8/3 

JnT JnV 

+ \Dn{p- Pb,p- Pb) - ( h{ex)m, 

^ JnV 

where £>„ is the Coulomb interaction 

Dn{f,g) = (/,5)ij-i(„r)' 
and the density p and the spin density m arc given by 

2,2 22 

p — ly^ + m = Vj^ — v_. 

The functional (5.1) is optimized under the normahzation constraint on the electron 
density 

(5.2) n-^ I p = / ul + v'i=Z. 

JnT JnT 

The Euler-Lagrange equations associated with (5.1) are given by 

(5.3) - Au^ + ^j.;/^ - ^4/^ + {V- h{ex))u+ = 0; 

(5.4) - Ai/_ + \v_!^ - ^i/!./^ + {V + h{sx))v- = 0; 

(5.5) - AV = An{vl+vl- pb), 

in nT with periodic boundary condition. Note that the normalization constraint 
(5.2) is contained in (5.5) as the solvability condition. The potential V is determined 
up to a constant by (5.2), the constant is fixed by (5.3) and (5.4) through the 
solvability condition of (5.2). 

For later use, let us also write down the Euler-Lagrange equations in £-scaling, 
which is just a rescaling of (5.3)-(5.5). 

(5.6) - e^A^ + ^«)V3 _ 1(^^)5/3 + (y. _ ^ 0. 

(5.7) - e^At/l + \{piy'^ - \{ylf'^ + {V^ + h)vl = 0; 

(5.8) -e^AV^ = An{{v%f + {vtf-pl), 

in r with periodic boundary condition. We have the scaling relations 

z/| {x) = u± {x/e) , {x) = V{x/e) , pi (x) = pb {x/s) . 

In analogy with the spirit of the Cauchy-Born rule for crystal lattices, we ex- 
pect that the electronic structure around a point xq to be approximately given 
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by the electronic structure of a crystal under constant applied potential with am- 
plitude h{exo)- As we have discussed in the last section, the electronic structure 
for the system with the constant applied potential is given by i'±,cb{-', h{exo)) and 
VcB{-]h{exo)). Therefore, the electronic structure constructed using the spirit of 
the Cauchy-Born rule is 

(5.9) 'y±{x) = iy±fiB{x;h{ex)), V{x) = VcB{x;h{ex)). 

One main result of this paper is that under the stability conditions, the electronic 
structure constructed by the Cauchy-Born rule gives a good approximation to a 
solution to the TFDW equation. In other words, one can find a solution to the 
TFDW equation that is close to the Cauchy-Born approximation. 

Theorem 2. Under Assumption A, there exist positive constants ho, Sq and 6, 
such that for any h G C°°(r), ||/i||i,oo(r) < ho and e < so, there exists a unique 
u = {1/+, v-jV) e {H^)^, with the properties 

• u is a solution to the Euler- Lagrange equation, 

F{u) = 0; 

• u is close to the approximation given by the Cauchy-Born rule 

\\u - ucb{x; h{ex))\\(H2)3 < 6e. 

We notice that once the solution to the Euler-Lagrange equation is determined 
as u = {i'±,V), the associated energy can be written as 



Jnr Jnr 

+ / V{p - Pb) - / h{ex)m. 

Jnr Jnr 

As a consequence of Theorem 2, the energy is well approximated by the Cauchy- 
Born rule, given by 

(5.10) InCB= [ EcB{h{ex})dx = ^ [ EcB{h{x))dx, 

Jnr ^ Jr 

with i^cB defined by 
(5.11) 

EcB{h)= ■ful% + u]_% + \Wu+,CB\' + \^t^-,CB\'- {- ^%B+^%B 

Jr Jnr 

+ j- Vcb(pcb - Pb) - hj-mcB, 

where i^±^cb = i'±,Cb(-;^) and similarly for other terms. Note that, at least for- 
mally, we have 



/ EcB{h{x))dx = mi( f 
Jr V^r 



(5.12) / EcB{h{x)) dx = inf / £;cB(m(a;)) - h{x)m{x) dx . 
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Here -Ecb(w) is given by the cell problem 



EcBim)= inf (-fuf^ + u]_' 





1 

2|f| 



Dt{p - Pb, P - Pb) \, 



where mtot = /j- v\ — v"^ is constrained to be equal to m. 

In terms of micromagnctics, the former term on the right hand side of (5.12) 
is the anistropic term of magnetization, and the latter term is the energy due to 
external magnetic field. Compared with the usual energy functional used in micro- 
magnetics [2,3], we do not have the stray field energy term (the nonlocal term) for 
the magnetostatic interaction and the exchange term. The reason that the nonlo- 
cal term is missing is due to the fact that the Thomas- Fermi-Dirac- von Weizsacker 
model only contains a local term of m. One can try to add a term that account for 
the magnetostatic interaction at the microscopic level, we would then obtain the 
nonlocal term in (5.12). However, this is only natural for a non-coUinear model, 
which will be studied in future publications. The reason that we do not have the 
exchange term in the energy functional is more fundamental. Since the scaling 
we consider only covers the smooth regime, there is no hope on the leading order 
to recover the exchange term which penalize change of magnetization on the scale 
comparable to the atomic length scale. One might hope to obtain the exchange 
term by a different scaling limit or going to the next order, for example, zooming in 
the region of domain wall. We would not go further in this direction in the current 
work. 



6.1. Matched asymptotics. In this section, we use two scale analysis to build a 
high-order approximate solution to the Euler-Lagrangc equation. 

It is more convenient to work with £-scaling in this section. The choice of scaling 
is also in agreement with standard homogenization problems. Let us recall the 
Euler-Lagrange equations under e-scaling. 



6. Two SCALE ANALYSIS 




-e^lW' = ^i:{{vlf + {vtf-pl). 



We take the following two-scale ansatz 



vL{x) = v-fi{x, xje) + EV-^\(x, xje) + e^v-^-2,(x, x/e); 
V^{x) = Vo{x,x/e) + eVi{x,x/e) + e'^V2{x,x/ e). 
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Substituting into the equations and matching orders, we obtain for the leading 
order 

(6.1) - A^iy+^o{x, z) + ^i^+(o(a;, z) - ^i^+(o(a;, z) 

+ {Vo{x,z) - h{x))iy+fi{x,z) = 0; 

(6.2) - A^zy_,o(a;, z) + ^vl/flix, z) - ^u^J^ix, z) 

+ {Vo{x, z) + h{x))y-fi{x, z) = 0; 

(6.3) - A^Fo(a;, z) = At:{vI q{x, z) + o(a;, z) - pb{z)). 

The solutions are given by 

(6.4) I'+fiix, z) = u+,cb{z; h{x)), f-flix, z) = z/-,cb(-2; h{x)), 
with the potential given by 

(6.5) Vo{x,z)=VcTi{z-Mx))- 
The next order equations are given by 

(6.6) - A^zz+^i - 2Vx ■ + - ^^l^li^+^i 

+ (Vo - h)i^+,i + Vii^+^o = 0; 

(6.7) - A^z/.,! - 2Vx • Vzi'-fi + ~ ^^-(o^^-.i 

+ (Vb + /i)z/-,i + Fi!/_,o = 0; 

(6.8) - A^Vi - 2Vx ■ V^Vb = 87r(!/+,oz/+,i + i^-,o^-,i)- 

Using the linearized operator Ch , we may rewrite the set of equations as 







(Ua\ 


(6.9) 




f-s 






\9iJ 



where 



By the regularity of j^±,o and Vq, it is easy to see that /±,i, gi G C°°(r,if^). 
Inverting £h by Proposition 4.1, we then obtain i^±,i, Vi, and 

i^±,i e c°°(r, Hf) Vi e C7°°(r, Hf). 



CAUCHY-BORN RULE FOR SPIN-POLARIZED TFDW MODEL 



17 



We also have the third order equations 



(6.10) -A^i/+,2-2V, 



.35 4/3 70 1/3 2 



3.11) 



- ^i'+fiiy+,2 - '^y+]i^yl,i + (Vb - h)y+,2 + + V2y+fi = 0; 

.35 4/3 70 1/3 2 

20 

^yZ^d^y^^i + {Vo + h)y-,2 + V^z^-.i + V2y-,o = 0; 



9 

^zy-,2 



2V, 



20 2/3 



(6.12) - A^Va - 2Vx • V^li - A^Vb = 47r(!/2 ^ ^ 2v+fiV+^2 + v'i^^ + 2i/_,oi^-,2). 



As the first order correction, we solve 



(6.13) 



^V+,2\ 




V~,2 


- f-J 




\92 J 



where 



70 



/+,2 = 2V^ • V^y+,1 + A^iy+fl - '-^y]lly\^^ + ^i^+'o'''^+,i - ^I'^+.i; 



20 



-1/3 „2 



27 

70 



1/3 2 

27"-'°"-'! 
1 



27 

^,,-l/3„2 
27" 



/_,2 = 2Vx • VzZ/_,l + Aa;ZV_^o 

52 = -^(2V^ • ^zVi + A^t/o) - o(''+ 1 + i)- 

OTT zi 

Therefore, (6.13) is solvable to give i^±,2 and V2 in C°°(r. Ff). 

This procedure can be carried on for even higher order terms. We omit the 
details here. 

Remark. There is an important difference between the Thomas-Fermi type of mod- 
els and the Kohn-Sham type of models considered in [5] . In the two scale analysis 
for the Kohn-Sham map developed in [5] , the macroscopic part of the potential on 
the leading order depends on the density on the order of e^, making the closure a bit 
unusual. Here, the macroscopic part (Vo) is determined on the leading order, and 
it imposes a constraint on the third order densities. In particular, as observed in 
[4] and [8] , the Coulomb potential in Thomas- Fermi type of models are determined 
locally to the leading order, while it is not the case for Kohn-Sham type of mod- 
els. This also leads to important differences in developing multiscale algorithms for 
these two type of models. 

6.2. Approximate solution. Let us take the approximate solution built in the 
last section: 

y+{x) = y+fi{ex,x) + sy+^i{ex,x) + e^y+^2{£x,x); 
^-{x) = y-fi{ex,x) + sy-^i{sx,x) + e'^y-^2i£x,x); 
V{x) = Vo{ex,x) + eVi{sx,x) + s'^V2{sx,x). 
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Here we have rescaled the functions into the units in which the lattice parameter 
is 1. 



\\T{y+,v.,V)\\(Li).<e\ 
{f+,f.,g)=T{v+,v.,V). 



Proposition 6.1. 

Proof. Denote 
We write 

= + f+a + /+,3 + f+A. 

where we have introduced the shorthand notation 

/+,2 = + + - ^i^+(o 

35 4/3/ , X 70 2 1/3 2 



27 



20 2/3. , ^ ,20 2 -1/3 2 



/+,4 = + e^V2V+,l + £'^V2U+,2. 

Here v and F are evaluated at (ex, a;). Since H/Hl^ < ll/lli~> it suffices to prove 
that < for i = 1,2,3,4. Since we have v+^2 e C°°{T,Hl), by 

Sobolcv embedding Hf C i°°(r), it is easy to see that ^ The 

estimates for f^ 2 /+.3 follow from Taylor expansion. Finally, the desired 
estimate of /+,4 is obtained from Sobolev embedding applied on Vi, V2, v+^i and 
y+,2- In summary, we have 

II/+IIl^<^'- 

The argument for /_ is completely the same as that for /+. Let us consider 

g = -e^^V^-A^({ylf + {yl_f-pl) 

(6.15) 

= 51+52+53, 
where we have introduced the shorthands 

51 = -e^A^^i - 2£3v, • V,V2 - e*A,y2; 

52 = 87re^z/+,ii/+,2 + 47re^i/^ 2; 

53 = di'Ke^V-^iV-^2 + 47r£^I/^ 2- 

It is clear that by analogous argument as above, we have the estimate ||5||l2 < e^. 

□ 
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7. Proof of Theorem 2 

Similar to [5,6], we use Newton- Raphson iteration to find a solution to the Euler- 
Lagrange equation in the neighborhood of the approximated solution constructed 
above. 

We will start with the approximate solution we constructed 
(7.1) u° = {uiy_,V'')GVn, 

where 

i^+ix) = v+fl{ex, x) + ez^+.i(£x, x) + e'^i^Ji.,2{sx, x); 

= v-^oi^x, x) + ei^__i(ex, x) + £^j^__2(ea;, x); 
V°{x) = Vo{ex,x) + £Vi{ex,x) + s'^V2{ex,x). 

Here the superscript is used to indicate the initial point for the Newton iteration. 
Note that we have rcscalcd the functions, so that ly'j. and arc defined in nT. 

We need several additional lemmas for the proof of Theorem 2. The following 
lemmas are proved under the assumptions in the statement of Theorem 2. 

Lemma 7.1. For any k > 0, there exists positive constants ho and Sq, such that 
for all s < eo and h with ||/i||ioc(p) < ho, we have 

\\CuO - £per||.Sf((L2)3) < K. 

Proof. By definition, given {lo+,uj-,W) G {H^)^, we have 









(5C+ 





(7.2) 


{^uO -^per) 


tA- 





5L 






[wj 




5v 












■ 




wj 







where the last equality serves as a definition and we have introduced the notations 

= f ((.« )V3 _ _ ^((,0)2/3 _ ^^yO_^_ y^^^^^ 

SC- = f ((.° )^/^ - <y - ^ ((.° )^/^ - .^%) + (y° + - Vpe.), 

and 

5u± = I/O - pg,. 

Note that 

diy+{x) = i'+^o{ex,x) + ei'+^i{sx,x) + s'^i'+^2{£x,x) — ly+^perix) 

= {'^+,CB{x;h{ex)) - zv+,cB(a;; 0)) + eiy+^i{ex,x) + £^z^+,2(ea;,x). 

By the smooth dependence on i^+,cb('; h) on the parameter h, we have 
\\i'+,CB{x;h{ex)) - iv+,cB(a;; O)||i,oo(„r) < ||/i||L~(r)- 
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The other two terms are of higher order in e, since :^+,i(ea;, x) and vj^^i{^x^x) are 
bounded uniformly in nF. Therefore, we obtain 

||'5i^+||L~(nr) < ||^||L~(r) +e. 

Obviously, the same estimates also hold for 8v^. Hence, 

llslU^ < ll'5i^+IU°°(nr)||w+||L2 + ||(5j^_||ioo(„r)||w_||L2 

<(|l/i|U»(r)+e)(||a;+|U. +||a;_|U.). 

The analysis for /+ and /_ are the same, let us study /+ . Using Taylor expansion, 
we have 



max u 

^e[l^+,mi„(a;),y+,„ax(2:)] 



1/3 



where iy+,niin{x) = min(i/!|(a;), i^+,per(a;)) and !^+,max(a;) = niax{u^{x),u+^per{x)). 
Since 

- J^+,per(a:)| < WS^lh^ < i|/i|Uo.(r) +e, 
we have for ho and eq sufficiently small, I'l (x) is bounded from above and also from 
below away from zero uniformly for x G nT. Hence maxyg[^^ min(x),i/+ max(x)] '^^^^ is 
bounded. Therefore, 

\iulf/^x) - {u+,p,,f/'{x)\ < \\hU^^r)+e. 



It follows that 



35 



(7.3) 



LI 



<(l|/^l|L~(r)+£)IKIlL^ 



Using similar arguments, we have 



(7.4) 



f(K)^/^-K,per)^/^)-+ 



< 



LI 



Compare the difference of and Vpor, we have 

V° - Vper = VcBix; h{sx)) - VcBix;0) + eVi{ex,x) + s^V2{sx,x). 
Analogous to the control of , we have 

||l^°-V^per|U~(„r)<||/l|Uoc(r)+£. 

Therefore, 

< (||/i|Uoo(r)+£)||w+|U2. 
We also have 

(7.6) \\5u+WUi < \\Su+Uoo^„r)\\WUi < (||/i|Uoo(r) + ^lU^ . 
Combining (7.3)-(7.6), we obtain 



(7.5) 
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Therefore, 

The conclusion of the Lemma follows. 



□ 



Corollary 7.2. There exists positive constants ho and eq, such that for all £ < £o 
and h with ||/i||L°°(r) ^ ho, we have 

ll^;;oi.Sf((i2)3,(H2)3) < 1. 

In particular, the bound is independent of n. 

Proof. Since 

we have 

~ ■^per(-^ + i^u" ~ ■^per)^per) 

if the right hand side is well-defined, where I is the identity operator. 

By Assumption A, there exists k > such that \\jC ||_5f((ilf2)3_(i2)3) < K 

implies 

IKO - £per)£pcr||^((//2)3) < 1/2, 

and hence T+ (£„o — Cpcr)i^pcr is invcrtiblc on {Hf-^f' with the norm of the inverse 
less than 2. Therefore, by Lemma 7.1, there exists ho and Eq sufficiently small that 

\\^Z^\\^{{Ll)'^AHir) ^ A^plr\\^((LlY,(Hl)3)- 

The corollary is proved by combining the above inequality with Assumption A. □ 
Lemma 7.3. If u,u' satisfy \\u — u^\\^^jjl)^ < 7^^) 11^' ~ — 7^^' ^^^'^ 

ll^u ~ ^u'\\j^((LlY) < C(7)e"^/^||u- m'||(//2)3. 
Proof. By definition, given e {I^iY ^ have 



(7.7) 







(5C+ 





5v^ 



















5u_ 














1 




■ 













where the last equality serves as a definition and we have introduced the notations 
35 / , , , ,4/o\ 20 , 



and 



First let us control g. By Sobolev inequality, we have 

ll^'^±||i,oo(„r) ^ ||^z^±||i?2(„r) =n^''^\\5v±\\Hi. 
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Hence, 

< ||(5zv+||i,oo(„r)||w+||i,2 + ||^!/_||ioo(„r)||w_||i,2 
<n3/2(||5i.+ || 

h^W'^+Wli + W^'^-Wh^W'^-Wli)- 

Similarly, we have 

\\6iy+W\\Li<n''^^\\6iy+\\H2jW\\Li. 
Now consider the term d£+u+ , we have 

\iu+r^\x) - iu'^r/'{x)\ < ^\u+ix) - u'^ix)\ max 
where z^+,min(a;) = min(iv+(a;), z/!^(a;)) and !^+,max(a;) = niax{i'+{x),i''_^{x)). Since 

we have for e sufficiently small, i^+{x) is bounded from above and also from be- 
low away from zero uniformly for x G nV. The same holds for and hence 
maxj,g[^^^,^(a.)_^^^^^(a.)] !/V3 is bounded. Therefore, 

It follows that 

The argument for the remaining two terms in 6jC+oj+ is similar, and we conclude 
that 

\\SC+co+Ui < C{j)n'^^u - u\hI)4^+\\li, 

and analogously 

In summary, we have ||(5£||^((l2 < C{'j)n^^'^\\u — w'||(//2)3 = C{'-f)e~^^-^\\u — 

U'||(H2)3. □ 

Lemma 7.4. If u,u' satisfy \\u — u°||(//2)3 < 7^^, \\u' — u^\\(h2)3 < -fe^ , then 

\\Cj{j^{u) - T{u') - £„o(« - u'))\\(Hir ^ Ci^y^Wu - n'||(H2)3. 
Proof. We write 

T{u) — F{u') — C^o(u — u') = I — £„o)(u — u') dt, 



where Ut = tu-\- (1 — t)u' . It is easy to see that for any t G [0, 1], we have 
\\ut - ■"°||(H2)3 < max(^||u - u°||(H2)3, \\u' - u°||(h2)3). 
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Therefore, using Lemma 7.3, we obtain 

\\J^{u) - J^{u') - Cuo{u - u')\\(^Ll)^ < sup \\Cut - Cy,«\\j^{{Lir-)\\u - u'\\(L2)3 

tG[0,l] 

<£-3/2C(7) sup ||«t-U°||(H2)3||w-u'||(z,2)3 

te[o,i] 

< C(7)e^/^||u-M'||(L2)3. 
We conclude using Corollary 7.2. □ 

Now we are ready to prove the main result. We will actually prove a stronger 
version of Theorem 2 with higher order error estimate. Theorem 2 is clearly a 
corollary of the following result. 

Theorem 2'. Under the same assumptions as Theorem 2, there exists a unique 
u = e {H^)^ such that 

• u is a solution to the Euler- Lagrange equation, 

F{u) = 0; 

• u is close to 

Proof. Consider the nonlinear iteration 

(7.8) u''+'=u''-Cjj^iu''), 

with initial condition uP given in (7.1). 

Consider the first iteration (fc = in (7.8)), we get 

ll"^ - "°II(H2)3 = ||/:^,,^J'(M")j|(^2)3 < ||£^o^j|^((^2)3 (^2)3)11 J"(u°)||(i2)3. 

By Proposition 6.1 and Corollary 7.2, there exists constant Ci such that 

llu^ -it°||(H2)3 < Cie^/2. 
Suppose we have proved for all k < ko, 

(7.9) ||u*=-u°||(^j2)3 <Cl^^ 
by the iteration scheme (7.8), we have for all k < ko, 

= -c-^{nu^) - nu'^-^) - £„o(«'= - u^-^)). 

Hence, by (7.9) and Lemma 7.4, 

For e sufficiently small so that C(Ci)£^/^ < 1/2, we then have 

k=0 
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Therefore, by induction, we have for all A: > 0, 

||^x'=-u°||(jj2)3 <Cl^^ 

and 

(7.10) \\u^+^-u\„.)^ < h^\u^-u^-\^.j,. 

Because of (7.10), the iteration converges to u* and 

T{u*) = 0. 

In addition, we have 

\\u*-u'\\^Hir<Cie^ 
The local uniqueness also follows easily from (7.10). □ 
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